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THE APPLICATION OF MATRIX METHODS TO COORDINATE 
TRANSFORMATIONS OCCURRING IN SYSTEMS STU DIES 
INVOLVING LARGE MOTIONS OF AIRCRAFT 
By Brian F. Doolin 


SUMMARY 


The purpose of this paper is to show the method and advantages of 
matrix algebra in setting up the geometric aspects of problems of airplane 
motion. Such aspects arise particularly when studies of systems which 
include aircraft are being made. The geometry is formulated by fixing 
quantities whose relative motions are to be studied, each in a coordinate 
system of its own. The various coordinate systems are related to 
other by orthogonal transformations in matrix form, end the parameters 
defining the transformations are found in terms of the dynamical variables 
of the problem with the help of the transformation matrices. The compact 
notation of matrix algebra permits a clear view of the geometry involved. 
Use of matrix algebra provides a routine procedure for computing the 
detailed expressions required in a particular problem. The first part 
of the paper discusses those aspects of matrix algebra required for use 
in orthogonal transformations. The second part shows how to use orthogonal 
transformations in matrix form by applying them in several examples. 


INTRODUCTION 


There are many problems currently "under study that concern the motion 
of one or more aircraft over relatively long periods of time. For example, 
long trajectories are involved in some studies of very high -altitude, high- 
speed aircraft. Another example is provided by current fire-control 
systems studies where the relative motion of two aircraft is involved. 

The elements of interest in such problems, in one case the position and 
attitude of the aircraft with respect to a nonflat earth, in another the 
relative positions and rates of the two aircraft, are referred to coordi- 
nate systems that undergo large changes in orientation. The formulation 
of these problems for study, say, on an electronic analog computer, 
requires the expression of these large orientation changes by means of 
orthogonal transformations. These transformations refer the coordinate 
systems, in which the elements of interest are imbedded, to an initial 
coordinate system whose orientation is fixed. 
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The equations that express the transformed components of a quantity- 
can he obtained in a simple fashion if the orthogonal transformations of 
the coor dina te systems involved are written in matrix form and the rules 
of matrix algebra followed. The matrices that perform orthogonal trans- 
formations are rotation matrices. A description of their properties and 
the algebra they satisfy - that is, the rules by which such operations 
as addition and multiplication are performed - can be found in many places 
(see, e.g., refs. 1, 2 , 3 , and 4). Their application to certain aspects 
of problems concerning aircraft is illustrated in references 1, and 8. 
One purpose of the present paper is to present the matrix method of per- 
forming orthogonal transformations in as simple a manner as possible, to 
describe its algebra, and to illustrate the algebra in use by several 
examples . 

The examples in this paper will serve a second purpose, namely to 
show that the complete geometry of the problem under study can be expressed 
in matrix form. By imbedding any element of interest in a coordinate 
system, and by following this coordinate system by means of matrices as 
it rotates with respect to all other elements of the problem which are 
fixed in rotating coordinates of their own, a single matrix expression 
for the over-all geometry is obtained. From such an expression the 
expression relating any two quantities of interest can be derived in a 
routine manner. Furthermore, since in any problem there are a number of 
ways of choosing an appropriate set of rotating coordinate systems, the 
compact notation of matrices which, provides a clear view of the geometry 
of the problem is especially valuable in permitting the examination of 
various possible choices of sets of coordinates, thus facilitating selec- 
tion of one set that may simplify the representation of the geometry on 
a computing machine. The matrix notation permits this ex ami nation to 
proceed without detailed computation having to be performed. 

Once the choice of coordinates has been made, the parameters that 
specify them, such as the angles cp, 6, and t that express the orientation 
of an airplane, can be determined in terms of quantities that occur in the 
problem solution, such as the airplane rotation rates p, q, and r. The 
examples in the present paper also show how the determination of these 
parameters, or their rates of change, can be obtained by matrix methods. 
Thus all the geometric expressions required in these problems can be 
obtained by the regular and routine application of matrix methods. 

A question of notation arises from problems such as those described 
in this paper. Not only are there several quantities of interest in each 
problem, but some of these quantities have to be expressed in several 
coordinate systems. The notation adopted in this paper to satisfy the 
requirements of this type of problem is described briefly in the text and 
more fully in the Appendix. 
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DISCUSSION OF BASIC HIOFERTIES OF ROTATION MATRICES 


The first part of the present paper is intended for those readers 
■who either are not at all familiar with the algebra of matrices or are 
not accustomed to thinking of matrices in terms of their properties of 
performing rotations. The fundamental matrix algebra, such as the con- 
ventions governing addition and multiplication, will be presented in a 
general manner. Orthogonal matrices (which perform only rotations, and 
therefore will be called rotation matrices) form a special group of 
matrices and have properties, not possessed by other matrices, which 
simplify their manipulation. Such properties will be pointed out and 
subsequently often used. 

A glance at the examples in later parts of this paper will show that 
the matrices required in a single problem may form a long chain performing 
rotations in succession. H3ie discussion of the basic properties of rota- 
tion matrices will begin with a consideration of matrices performing a 
single rotation. Then the properties of matrices of successive rotations 
will be discussed. 


Single Rotations 


The rotation matrix performs an orthogonal transformation on some 
quantity. This quantity may be a vector or a set of vectors that define 
a particular coordinate system. The coordinate systems considered are 
three-dimensional and Cartesian. Three vectors are required to define 
such a coordinate system. In this paper, a triad of vectors defining a 
coordinate system (or set of coordinates) are of unit length, at right 
angles to each other (mutually orthogonal), and oriented with respect to 
each other in such a way as to form a "right-handed" system. 


Rotation matrices perform the operation of rotation on these quanti- 
ties i that is, the quantity obtained by transformation differs only in 
orientation from the original quantity. Not only does It retain its size, 
but, if a coordinate system is the rotated quantity, the unit vectors of 
the new coordinate system are mutually orthogonal in the same way as were 
the unit vectors of the old. 


Special forms of rotation 
matrices , - The matrices of simplest 
form are those that use one or 
another of the three unit vectors 
as an axis of rotation. In the 
sketch adjacent, a new coordinate 
system is obtained from the old by 
a rotation_through the angle 9 
about the 1 direction. Several 



'o ' 'n 
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items c&n be noted in the sketch. The three uni t, vectors are labeled 
1, 2, 3. The original^set has the subscript o; the new, the subscript n. 
The rotation with the 1 vector as an axis is said to have a positive sense 
if it has the sense_of rotation of a “right-hand screw" which is advancing 
along the positive 1 direction. 

If an old and a new set of coordinates are related as shown in 
sketch (a), the components of a vector (or the coordinates of a point) 
in the new set are related to its components in the old by the following 
set of equations: 


xin = xio 


X 2 n = x^cos cp+x^ 0 sin cp l 

xa n = -Xaosin cp+X 3 0 c ° s cp 


(la) 


which is the same relation obtaining between the new unit vectors and the 
old 


In = lo 

2n = 2 0 cos cp+3 0 sin cp 
3 q = -2oSin cp+3 0 cos cp 


(lb) 




¥ 


This system can be put into the following tabular form from which the shape 
of the matrix emerges : 


\xio 

x 20 

*30 


x m 

1 

0 

0 


x 2 n 

0 

cos cp 

sin cp 

( 2 a) 

Xsn 

0 

-sin cp 

cos cp 



How to get equations (l) from the table is clear: one reads along a row, 

multiplying each member of a row by its column heading. The arrangement 
can be changed to read like an equation if proper conventions are adopted 
to multiply the elements of the table by the correct quantity to give 
equations (l): 


xin 

= 

0 

0 


x 10 

x 2n 

= 

0 cos cp sin cp 


x 20 


= 

0 -sin cp cos cp 




(2b) 
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With the proper conventions adopted, this table is an equation ■which can 
he symbolized thus : 

[new components] = (9) [old components] (2c) 

where (9) is the matrix with elements: 

“10 0 “ 

0 cos 9 sin 9 
_0 -sin 9 cos 9 


Equality indicated between two matrices holds if and only if 
corresponding elements are equal. Corresponding elements of two matrices 
are the elements in the same row and column in the two matrices. The 
convention used in multiplying matrices can be determined by comparing 
equations (la) and ( 2 b). In equation ( 2 b) there are three matrices 
whose symbols are [new components], (9), and [old components]. The first 
and third, column matrices, have three rows of elements all forming a 
single column. The (9) matrix has nine elements grouped into three rows 
and three columns. Equations (l) indicate that the product ( 2 b) or ( 2 c) 
stands for three equations. To obtain the first of equations (l), the 
elements of a given row of (9) are multiplied by each one of the element s 
of the following column matrix, and the products added. For example, 
for the first row in (9), the element of column 1 is multiplied by the 
element of row 1 in [old components]; the element of column 2 in (9) is 
multiplied by the element of row 2 in [old components]; element of 
column 3 in (9) is multiplied by the element of row 3 in [old components]. 
The sum of these products is equal to the element of the first row in 
[new components]. The procedure is repeated for each row of (9) until 
all three rows of [new components] are obtained. 

The three simplest rotation matrices together with a sketch of the 
operation each performs are shown below: 
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3 0 =3 n 


Sketch (c) 



The three special forms of rotation matrices given next to the sketches 
are the operators most commonly employed because they are simple in form 
and easy to visualize. It is not necessary for one to know their exact 
forms offhand, since they are easily derived. 


General form of rotation matrix .- Each of the three rotation matrices 
in equation ( 3 ) is a special case of the general form of rotation matrix. 
Ohe derivation of this general form follows easily from certain consider- 
ations concerning vectors. (The notation used below and in subsequent 
sections of this paper is explained in the Appendix.) A vector is a 
geometric or dynamic quantity, containing a magnitude and a sense, which 
exists independently of any description of it in a particular coordinate 
system. Thus, if v is some vector, say, a velocity, then 

V = V 1 n l^+V pn 2 3 = ViqI +V2q2 +V 3 o3 
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■where the subscript and superscript n’s and o’ s refer to new and old 
coordinates, respectively (the new having “been obtained from the old 
coordinates by a rotation), and the subscript numbers indicate the com- 
ponent of the vector. Hie dot product of V with 1,2, and 3 turn 
gives 

j— T l — O — £1 — “O T7"H ^rO - H 

*V*1 = Vyn = Viol *1 +Vso2 *1 +Vqo3 *1 

V*2 n = V2H = Viol 0 '^ +V20? 0 * 2\v 30 3°-~£ a 

v-I n = v 3EL = v 10 r 0 *3 n +v 2 o2 0 ‘3 n +'V'3oT 0 *3 3:1 


— n — n 

since 1 *1 


-n — n — n — n -n — n — n — n — n — n 

2 *2 = 3 '3 =1 and 1-2 =1*3 = 2 *3 = 0. 


Hiis set 


of equations can be represented in the matrix form 


— — 


_o _n _o _n _jo _n 


— 

Vm, 


1 *1 2 *1 3 *1 

— 0 — n — 0 — n — 0 — n 


Vic 

v 2n 

s 

1-2 2 *2 3-2 

— 0 — n — 0 — n — 0 — n 


V 2 c 

v 3 n 


1-3 2 .3 3-3 


V 3C 


w 


Hie 3 x 3 matrix in ( 4 ) is the desired general expression. Since the vectors 
in the matrix are of unit magnitude, the elements in the matrix are simply 
cosines of angles. Each of the three rotation matrices in (3) is a special 
case of that in ( 4 ) obtained by making one of the unit vectors of the new 
coordinate set coincide with one of those of the old. Hie sines of angles 
in (3) are the cosines of 90° plus or minus those angles. 


Successive Rotations 


Matrix multiplication, - In writing out the terms resultin g from the 
multiplication of a three row, three col umn matrix (called a 3x3 matr ix ) 
by another 3x3 matrix, one follows the same procedure as that described 
in the previous section. Hie product is a 3x3 matrix, rather than the 3 x 1 
matrix obtained above. Hie procedure described furnishes the elements of 
one column of idle product. For example, let (C) = (A)(B): 
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Cxi 

c 12 

C 1 S 


a ll 

a 12 

a lS 


in 

bi2 

bis* 

C 2 1 

c 22 

C23 

= 

a 21 

a 22 

a 23 


b 2 i 

b 22 

b 23 

C31 

c 32 

C33 


asi 

a 32 

a 33 


bsi 

b 32 

b 33 _ 


[ 5 1 ib 1 i+a 12b 2 i+a 13b 3 1 
a El^ ll+ a 22^21+a2S^3 1 
a 3 lt> 11+032^2 1+333^3 1 


a Ilk 12+ a 12^ ) 22+ a 13^32 

a 2 lb 1 p+appb 22 +a 23 b 32 
®3 lb 12+ a 32b22+ a 33bs2 


a llb 13 + a 12^23+3 13^33 
a 2lb l3+ a 22b23+ a 23"b33 
®3 lb l3+ a 32b23+a3QTD 3 3 


The subscripts on the various elements indicate the position of the ele- 
ments in the matrix. The first subscript indicates the row, the second, 
the column the element occupies. Since the elements of two matrices are 
equal if the matrices are equal, nine equations are contained in 
(C) = (A)(B). Two of the nine are: 


cn = a ! ib ! 3+a x^b 2 i+a 13! 3 1 ; c 12 = a x ib 12+a 12 b 22+a i3b 32 

3 

All nine equations are contained in the expression c-j^ * 

3=1 

This expression also describes the method of multiplication, for, to 
obtain the element in a given row (i = 1, 2, or 3) and column (k = 1, 2, 
or 3) i& (C), each of the j column elements (j = 1, 2, and 3 in turn) 
of row i in (A) are multiplied by each corresponding row element j 
(the same j =1; 2, and 3 in turn) of column k in (b) and the resulting 
terms are added. 

This multiplication procedure is general. 1 It therefore applies to 

two special, multipliers, the unit matrix and a constant.' The unit matrix 

corresponds to the number 1 in arithmetic which has the property that the 

product of the multiplication of any number by 1 is identical to the orlgi- 

nal number. This property can be stated operationally in the following 

iThe multiplication procedure is general provided the matrices are 

conformable (ref. 1, p.6). Conformability means that in the product 

3 

(C) = (A)(B) or cjjs =^ a ijb ( jk there must be as many elements with 

3=i 

index j in bj^ as there are in ai j . 
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way: Unity is that multiplier which operates on a multiplicand in such 
a fashion that the resulting product is the original multiplicand. The 
matrix with this property is the uni t, matrix, which has the form: 


10 0 
0 10 
0 0 1 


It can be written symbolically as (i), and is the "diagonal matrix" all 
of whose main-diagonal elements are unity. 

The product of a constant times a matrix is a matrix each element 
of which is the constant times the corresponding element of the multipli- 
cand matrix. Thus, If (C) = k(B), then ay = kbj_ ^ . The equation can 
also be written 


t>iibi 2 bi 3 


k 0 6 " 


tub 12b 13 


= 

0 k 0 


t 2 lb 22^23 

^ 3 xh 32 t 33 _ 


0 0 k 


t 3 lb 3 2 b 3 3 _ 


A constant can therefore be called a diagonal matrix with equal and 
constant elements. 

Matrix addition. - The sum of two matrices is defined as the matrix 
whose elements are the sums of the corresponding elements of the added 
matrices: let 



a XX a X 2 a X 3 


txxt>x 2 tx 3 

(A) = 

a 2 X a 22 a 23 

and (B ) = 

1 2 X^ 22^23 


a S X a 32 a 33 _ 


] 23 xb 32 b 3 3 _ 


Then (A) + (B) is the matrix: 

B-ii+bn a i g+b i g ais+bi3 

a 21+b2IL apg+bgp &23+b23 
_^ 3 1+^31 a 32 +t >32 S- 33 +b 33 

This identity can be written in the form (A)+(B) = (A+B), 
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Geometric Interpretation of rotation matrices, - Hie rotation matrix 
not only partakes of all the properties of matrices given above, hut also 
may he interpreted geometrically as an operator that rotates a coordinate 
system through 6 ome angle about some axis given in the coordinate system. 
Multiplication of two such operators corresponds to two successive rota- 
tions. If the operation depicted in sketch (a) is followed by another 
identical operation, the total operation must consist in a similar rotation 
of twice the magnitude of either one. That the algebra agrees with the 
geometry is seen by multiplying the rotation matrices together: 



1 0 0 ■ 


1 0 0 ~ 


0 cos q>e sin tp 2 


0 cos epi sin tpi 


0 -sin cp 2 cos ip 2 


0 -sin tpi cos 



1 

0 

0 



0 

0 

m 

0 

cos qj 2 cos cpi-sin cp a sin cpi 

cos cpjjsln q)i+sin cp 2 cos tpi 

- 

0 

COBflfe-KPi) 

sinCqjjj-KPi) 


_o 

-sin qigcos qpj-cos <p 2 sin <p x 

-sin (p 2 sin <Pi+eos <p 2 cos (p,. 


0 

-sinCifewpi.) 

cos((p 2+ <Pi) 


The order of the matrices Implies that the operation (cpx) is first 
to be performed on whatever is to its right, followed by operating by (92) 
on everything to its right. Hie order in which two matrices are written 
Is usually important. (When a number of finite rotations are performed 
about the same fixed axis, the order is not significant.) The following 
examples of rotation about two different axes illustrate that the order 
in which the rotations are performed is important. 


1 

0 

0 ~ 


"cos 0 

0 

-sin 

9 

0 

cos 9 

sin 9 


0 

1 

0 


0 

-sin 9 

cos 9 

— 


sin 0 

L. 

0 

cos 

9 


"cos 9 

0 

-sin 


1 

0 

0 “ 

0 

1 

0 


0 

cos 9 

sin 9 

sin 9 

0 

cos 0 


0 

-sin 9 

cos 9 


r cos 

9 


0 


-sin 

9 

— 

sin 

9 

sin 9 


cos 

9 

sin 9 

cos 

0 

cos 

9 

sin 0 

-sin 

9 

cos 9 

cos 

0 

cos 

e 

sin 

0 

sin 

9 

-sin 0 

cos 

9 

0 


cos 

9 


sin 9 


sin 

e 

-cos 

0 

sin 

9 

cos 9 

cos 

9 
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If the cp’s and 0 's are respectively equal in the two products, then 
the product matrices are not equal. 

Matrix inversion, - After a new coordinate system has "been obtained 
by rotation from an old, it is clearly possible to regain the old by 
another rotation. The second rotation is said to be the inverse rotation. 
If the matrix of the first rotation is (cp ) , that of the inverse is denoted 
by (cp)* 1 . In symbols, if 

[new components] = (9) [old components] 

then 

[old components] = (cp) -1 [new components] 

It follows, by multiplying the first equation by (cp ) -1 and the second by 
(cp) that (q>) (qp ) “ 1 = (l) = (cp) - 1 (cp) where (i) is the unit matrix. Written 
out, the matrix (cp) and its inverse have the following form: 



1 

0 

0 


1 

0 

-*> 

0 

(9) = 

0 

cos 9 

sin 9 

11 

H 

1 

*&■ 
' — * 

0 

cos 9 

-sin 9 


0 

-sin 9 

cos <£ 


0 

sin 9 

cos 9_ 


Upon examining the two matrices, one sees that the inverse matrix can be 
obtained from the (cp) matrix by exchanging the angle -cp for every cp. 
This property of rotation matrices, which is not possessed by other 
matrices, again confirms that the results of the matrix algebra agree 
with the geometry, for the inverse rotation geometrically is the same as 
a rotation in a negative sense. 

Inversion of a rotation matrix can also be achieved by transposing 
its elements, where the transposed matrix is obtained by exchanging ele- 
ments across the main diagonal of the original matrix; that is, if (B) 
is the transpose of (A), then bpj = aji, for example, b 12 = a 21 . In 
(cp), since this operation replaces sin cp by -sin cp and vice versa, it 
is equivalent to replacing cp by - 9 . 

Reversal property. - Suppose there has been a sequence of rotations, 
for e-ramp l e, a. (cp) followed by a (©) rotation. A consideration of the 
geometry indicates that to return to the old coordinates from the new, 
one must rotate first through (- 0 ) or (0)" 1 , then through (-cp) or (cp)” 1 . 
It is easy to show that again the algebra and the geometry agree. First 
of all the product operation ( 0) (cp) is itself an orthogonal rotation 
operator. Hence [(0)(cp)] -;L [(0)(cp)] = (i). But since ( 0 )“ 1 ( 0 ) = (i) and 
(cp)" 1 (cp) = (i), it must be that (cp) - 1 (0) _ 1 (0)(cp) = (i). As a result, 
(cp)” 1 ^) -1 = [( 0 )(cp) ] _1 . In words, the inverse of a product of matrices 
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is the product of the inverses of the matrices taken in reversed order. 
It is worth repeating that the property (-cp) = (cp) -1 is not a general 
property of matrices, hut belongs to those matrices that perform only a 
rotation. 


EXAMPLE APPLICATIONS 

Rotation of a Vector About an Arbitrary Axis 


For a first example, the form of single rotation matrices will he 
generalized. The rotation will no longer he performed about one of the 
old coor dina te directions, hut about an arbitrary axis, that is about an 
axis n of uni t length with components ni, n 2 , and n 3 in the given 
coordinate system. Another difference from the previous discussion is 
that in this example a vector p is to be rotated about the axis n 
through some arbitrary angle cp into a new position q. Heretofore, the 
coordinate system was rotated, and to transform a vector meant to express 
the same vector in two coordinate systems with the same origin but with 
one system rotated with respect to the other. This difference, however, 
reduces to a matter of sign, for to rotate a coordinate system through a 
positive angle is entirely equivalent to rotating a vector in the (now 
fixed) coordinate system through the same angle but in the opposite sense. 

The quantities given in the problem are illus- 
trated in sketch (d) . A unit vector n has the 
components m, n 2 , and 103 in the coordinate system 
formed by the triad of unit vectors 1, 2, and 3 . 
Therefore n = nil+n 2 2+nQ3 and ni 2 +n 2 2 +i% 2 =1. It 
is desired to rotate another vector JJ (with given 
components p x , p 2 , and P 3 ) through some given angle 
cp into the new position shown in the sketch as the 
vector q. The problem is to find the new components 
q x , q 2 , and <33 in terms of cp and the components of 
p and n. 

This problem has been solved both by matrix 
methods (ref. 1, p. 253 ; and ref. 3> P. 9^) and by 
pure vector methods (ref. 2, p. 168 ). The method here will use only the 
matrices already discussed, and will illustrate the approach adopted in 
succeeding examples in this paper. The general procedure is as followsj 
The vector n will be thought of as one of a triad of unit vectors forming 
a coordinate system. The vector £ will be expressed in this coordinate 
system which will then be rotated in a negative sense through the desired 
angle cp. The resulting new components of p will then be returned to 
the fixed coordinate system by an inverse transformation. These will be 
the desired components of q. 
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Let [(sKt)]” 1 be the matrix transforming p from the fixed into 
the n coordinate system. Then 


where 


q.i 


"pi 

12 


P2 

5s_ 


i>8. 


(5) 


(fl)(t) = 


cos 0 0 -sin 0] 

0 10 
sin 0 0 cos 0 


cos ijr sin ijr 0 
-sin ijr cos ijr 0 
0 0 1 


and 


(q>) = 


cos 0 cos ijr cos 0 sin ijr -sin 0 
-sin ijr cos ijr 0 

[sin 0 cos ijr sin 0 sin i|r cos 0] 


10 0 
0 cos cp sin cp 

to -sin cp cos cp 


The angles 1^ and 0 can "be eliminated hy means of the components of n 
in the fixed coordinate system. Since 


nx 


"cos 0 cos ijr cos 0 sin ij- -sin 0 


1 

n 2 

= [ (©) (ijr) ]n = 

-sin ijr cos ijr 0 


0 

n 3 _ 


sin 0 cos ijr sin 0 sin ijr cos 0_ 


_ 0 _ 


nx = cos 0 cos ijr; n 2 = -sin ijr; 113 = sin 0 cos ijr. From the latter 
expressions one gets: cos ijr = n^+ne 2 ; cos 0 = nx/\(ni 2 +n3 2 ; 

sin 0 = nQ/Vni^+nQ 2 . Therefore, 
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and 


[(e)(t)] - 


ni 

-nxn a 

-“a -] 

n/ nx 2 +h3 £ 

*7 ni 2 +nQ 2 

n 2 

ni 2 -^ 2 

0 


-n s n 3 

HX 

\/nx 2 +na 2 

n/ nxS+ns^ 


C ( 0)(^)] _1 


! ni 

s/rii £ +n3 2 

-ns 

*ini 2 +n q 2 



Performing the operations indicated in equation (5) one obtains a matrix 
with the following elements: 

The elements of the first column 


nx 2 + 


ni&tig^cos cp ns 2 cos cp 
(ni e +nQ a ) (m 2 +n^) 


= | ni 2 ng 2 cos q>fn 3 2 cos cp+n 1 s(n 1 2 +n a 2 )cos £ . ni a cos qj 

(nx 2 +na 2 ) 


= nxnx(l-cos cp)+cos 9 
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n i3-L2-D-in 2 cos cp +nssin 9 

ELiria + nin3n g Seos cp-niP3Cos cp _ n g ns 2 sin cp+n g nx 2 sin (p 
ni 2 + n 3 2 tlx^+^s 2 


= niPa+ £ 
n x 2 +n3 2 


-n g sin cp 


The elements of the second column 


= nins(l-cos cp) - n g sin cp 


nin g (l-cos cp)-n3sln cp 

n g 2 +(nx 2 +n a 2)cos cp = n^fl-cos cp)+cos cp 
n g na(l-cos cp)+nisin cp 

The elements of the third column 


nxiiQ + n i n Q rL g gcos T-^iriacos cp + n g na 2 sin cp+n g n a 2 sin cp 
niS+na 2 n-i^+ngS 


= nxnoCl-cos cp)+n g sin cp 


n g n 3 (l-cos cp)-nxsin cp 

1132+ D g g n 3 2 coB cp+m 2 cos cp 
nx 2 +n 3 2 


= n32(l-cos t) i n 2 gn 3 2cos qH-nx 2 cos cp+na 2 (nx 2 +r>Q g )cos cp 

Qi g +n 3 2 


= n 3 n 3 (l-cos cp)+cos cp 
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With these results, equation ( 5 ) can he written in the form 




"11^1(1^006 cp)+cos 9 mn 2 (l-cos cpJ-nQsin cp nirts(l-cos cp)+n 2 sin cp" 


Pi 

qg 

= 

nin 2 (l-cos cp)+n3sin cp n 2 n 2 (l-cos cp)+cos cp_ ^^(l-cos 9)-ni8in cp 


P 2 

%_ 


niiigCl-cos qj)-n 2 6in 9 r^n^l-cos cp)+nisin cp rigi^l-coB 9)4-006 9 


53 - 



nini n x n 2 nix^" 


”0 -rig n 2 



~P2 

(l-cos 9) 

n 2 ni n 2 n 2 n 2 n3 

4 -cos 9 (l) 4 -(sin 9) 

Hq 0 -n x 



P2 

> 

_nani n 3 n 2 nan^ 


-n 2 nx 0 





where (i) is the unit matrix. To see that expressions ( 3 ) are contained 
in equation (6) one must recall that whereas equation (6) expresses the 
rotation of a vector in a fixed coordinate system, expressions ( 3 ) rotate 
the coordinate systdm with the vector fixed. Since the two operations 
are opposite in sense, the angles in the matrices will he opposite in 
sign. With this in mind, after letting S' equal first n x l, then n 2 5, 
then na^ and letting cp equal first -9, then -9, then -i|r, one ohtainB 
the expressions ( 3 ). 

A valuable hy-product of equation ( 6 ) is obtained from considering 
the vector q - p generated during a small time interval by a rotation 
through the small angle 9. If, in the limit as At-»0, 9 /At-Xp = w, 
then 


lim hl£2SL£l « 0 aid 
At->o At 


lim = 

At->0 At 


P 


where p is the velocity of the point p as the vector p rotates about 
n with angular speed m. Equation (6) becomes 


llm 

At-=*-o At 



nini-1 

njiia Q-i^s 


pi 



0 

“"3 

n 2 


Pi 

Um Lk 

Ab- 4 »-o At 

nsD.x 

H2n2-1 narig 


Pe 

+ lim 
At-»o 

sin cp 
At 

ns 

0 

-Hi 


P2 


DqUi 

L_ 

ngii2 nan® 


Ps 





_n 2 

ni 

0 


* 
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:g 


or 



where cji, and U& are the components of the angular velocity vector 
to and u X p is the vector or cross product of the vectors cj and p. 


Location of an Airplane on the Surface of the Earth 


One problem of some current interest concerns the study of the motion 
of an airplane over a significantly large sector of the surface of the 
earth. With a simplifying assumption, the analytical expressions of the 
geometry required in this problem provide a good example of the use of 
matrices in orthogonal transformations. The reason for the need of 
geometry here is that the airplane equations of motion require specifi- 
cation of the gravity vector direction which changes as the airplane moves 
over the earth. For the example, the gravity vector is assumed to be 
directed toward the origin of a set of coordinates which is located at 
the center of the earth (called earth coordinates). 

The expression of the direction of the gravity vector with respect 
to the airplane requires knowledge of the position and orientation (that 
is, heading) of the airplane with respect to the earth coordinates. This 
knowledge is determined from the time history of airplane velocity, and 
the specifications of the airplane’s initial position and velocity. The 
geometry is set up by establishing a set of earth coordinates which 
rotate in such a way that the moving airplane is located on one of the 
coordinate directions. In the sketches below, which are intended to 



Sketch (e) Sketch (f) 
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define the parameters used in the geometry, earth coordinates are labeled 
e, and airplane body coordinates are labeled a. 

Sketch (e) shows how "present earth coordinates" are derived from an 
original set of earth coordinates: the sequence of rotations, each of 

which is made in a positive sense, is first through the angle Y, then 
through the angle 0, finally through the angle 4. Each rotation has the 
matrix associated with it which is given by the expression with the cor- 
responding lower case Greek letter in expressions ( 3 ). Three similar 
rotations through angles \[r, 0, and <p serve to determine the present 
airplane coordinate orientation with respect to the airplane' s original 
coordinate system. Sketch (f ) shows the relation between the original 
earth and original airplane coordinate systems . The original airplane 
coordinates are derived from the original earth coordinates by a 0-like 
rotation about the 2 e ° axis through 90° in a negative sense. Two other 
angles must be specified, which relate the airplane's velocity vector to 
airplane coordinates. These angles are shown in sketch (g). 



Sketch (g) 


expressed V in earth coordinates 
and 6. 


The sketch shows that the airplane 
coordinate system is obtained from the 
airplane velocity coordinate system 
(airplane wind coordinates) first by 
a positive rotation through the angle 
-P (p is the angle of sideslip), then 
by a positive rotation through the 
angle of attack a. These angular 
designations are conventional. 

With the geometry so_set up, the 
gravity vector is g = -gl_; the air- 
plane velocity vector is V = vl V j and 
the airplane position vector R = Rl. 
There remains the need to express g 
in aircraft coordinates and, having 
, to determine the angular rates Y, 0, 


The equations connecting the various coordinate systems can be 
expressed as follows: 


[present airplane coordinates] 


[present airplane coordinates] 
[present earth coordinates] 


(cp)(@)(>|f)[ initial airplane coordinates] 
(q))( 0 )(i)f)( - 90 °) [initial earth coordinates] 
(a) (-p) [present wind coordinates] 

( 4 )(8)(Y) [initial earth coordinates] 
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From these equations one has 


[present airplane coordinates] 


= (<P)(0)(t)( - 9°°H($) (9)(?)3“ 1 [ present earth coordinates] 


( 8 ) 


or 


and 


Sa - (cp)(e)(t)(-90°)[(®)(8)(?)r 1 g e 


[present earth coordinates] 


= j^(9)(0)(lO(-9O o )[ (&)( 0 )(¥)]“ 1 "j- (a) (-£)[ present wind coordinates^ 


>(9) 


or 


Ve =■ {(<P) (9) (t) ( -90°) [ (» ) (9 ) « ] -^""'Ca) 


where the subscripts on the vectors imply that the components of the 
vectors are written in the coordinate system indicated by the subscript. 
The superscript -1 on the product matrix of equation ( 8 ) in brackets 
means that the inverse of the product [ ($>) (0) ('£')] is desired. It will be 
recalled that [ (4>) (0) (y ) ]” x = = (-¥)(-©) (-«) . Since the 

whole product matrix of equation ( 8 ) is used in inverse form in equa- 
tion (9), it will be convenient to obtain it or its inverse as a unit 
first and then form the product with (a) (- 3 ). 

Equations ( 8 ) and ( 9 ) can be simplified by setting the angle <D = 0, 
or the matrix ($) = (i). Only two parameters are necessary to specify 
the orientation of a vector in each of the two equations. Again because 
only two parameters are needed, (a) and ( -p } are sufficient to express 
the airplane velocity vector in airplane coordinates . 

Equation ( 9 ) is: 


^ie 


“c© 

0 

-E0- 


~cY 

sT 

<r 


"0 

0 

-I 

v ae 

- 

0 

1 

0 


-sY 

cY 

0 


0 

1 

0 

Jse 


6© 

0 

=6 


0 

0 

1 


1 

0 

0 


ct -st <f 


cS 0 si" 


0 

0 

— 1 

sy cy 0 


0 10 


0 cq> -sep 

_0 0 1 


0 C8_ 


_0 S<p ©£ 



ca 


0 


-sg 0] 
eg 0 
0 1 
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where c stands for cosine and s stands for sine, and where 


(-90°) 


,o\-i 


“c(-90°) 

0 

s(- 9 o°r 


0 

0 -r, 

0 

1 

0 

= 

0 

1 

0 

-s(-90°) 

0 

c(- 90 °)_ 


_! 

0 

0 _ 


vie 

V2e 

v se. 


‘^s'0c0ci|r 

-s0( -si|fcq>+scps0ci)f) 

-se(scpst+c<ps0cnjr) 

+c0s¥c0s^ 

+c0s ¥ ( cqpcijf +s cps 0s\|r ) 

+c0s?( -scpcnjr+c£ps0silf ) 

+c0c¥s 0 

-c0c¥scpc0 

-»c0c'5'«pc0 

clcdBif 

e’3f( c<pci)f+scps0sijf ) 

cHf( -scpc++ccps03^r) 

-bUb e 

+s¥scpc0 

+s¥(ccpc0) 

cQcQcty 

c©( -styccp+scp80c\jr) 

ce ( aips^+ccpB 0 cnjf ) 

4 -s 8 sYc 0 si|r 

+s©s¥( cq)c\Jf+scp80si)f ) 

+80s¥( -scpe\|r+cq>60si|r) 

+B0c¥s0 

-s©cYs(pc0 

-S0c¥ccpc0 


cacp -casp -so | [v 
sp cp 0 0 

sacp -spsa ca 0 


(10) 


Because of the presence of only one el em e n t in the column matrix on the 
right, only the first column of the product of the two 3x3 matrices is 
required. Performing the multiplication, one gets: 

= cacp [ -s0c0c^r+c©s1 r s\|fce+c©c , ifse]+sp [ -s0( -s\|rccp+scps0ct )+c06^( ccpcijr+scpsest) - 
c0c'?scpc0]+sac!3[-s0(scpst+ccps0cijf)+c©sY(-scpCT|r+ccps0st)-c0cfccpc0] 


s cacp [ cYcOstJt-s Y s 0 ]+sp [ c¥( ccpcTjr+scps0ST|;)+s Yscpc0 ]+sac3[ c¥( -scpajf+ccpsSsiJr ) + 
v 

sfccpcO] 

= cacp [ c©c0ct+s©BYc0st+B©c'Fs03+sg [ c©( -STjfccp+scps0ci{f )+s0s¥( c<pc\if+s9s0st ) * 
v 

s©cYs<pc0]+sac3[c©(scpsi|f+ccps0ct)+^sl(-scpci)f+ccps0STjr)-s0cIcq)c0] 
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To express the gravity vector in airplane coordinates , the inverse of the 
first (3x3) matrix on the right side of equation (10) is waited. This 
inversion is performed hy transposing the matrix. Since g = -gl 6 , only 
the elements of the first row of the matrix in equation (10) are needed. 
The resulting equations are 


gxa = (-g) [ -S0c0c\jr+c0s¥c0sf+c0c¥s0] 

g 2 a = (~g)[ -s©(-silfcqM-scpsSci(f)+c©sl?(c9a|f+scpsesTjf)-c0c?s<pc03 
gsa = (“g)[-s0(sq)s'>lf+c9s0c^)+c©s1f(-scpci(r+ccps0s^)-c0c'fccpc0] 


The rotation parameters cp, 0, i|r, 0, I, which are continuous functions 
of time, remain to he determined. It is assumed that the rate of rotation 
of the airplane, fi a = pl^q^+rj 8, , is know continuously from airplane 
equations of motion. The letters p, q, r are the conventional terms 
for the components of fi a in airplane coordinates. The rates q>, §, ■£ 
will he expressed in terms of p, q, r. In the same way, 0 -will 

he expressed in terms of the components of angular rate about the earth 
axes, fif e , fif e , and fife. In turn, the angular rate of rotation of fi e 
is obtained from the knowledge of the airplane velocities in earth coordi- 
nates which are subject to the two conditions that R = R1 , that is, that 
the coordinates so rotate that the coordinate vector l e remains pointing 
at the airplane; and that <J> = 4 =0. 

One convenient method for obtaining the expressions for cp, &, and i]r. 
is given below. Other methods can be found in references 2 and 4. The 
angular velocity vector can be expressed in the forms 

fi a = p^ a +q2 a +r3 a =cpI a + §2 a, +« ao, (ll) 


where the coordinate systems indicated by primes are intermediate 
coordinates defined by 


~I a " 


-_a»- 


”“ a o r 


~l&o~ 

2 a 

-It 

= ( 9 ) 

2 S ’ 

—a* 

l_3 J 

= (cp) (©) 

2 a o' 

_3 ao - 

= (<p)(e)(t) 

2 a o 

■jB-O 


( 12 ) 
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Because of the form of the various rotation matrices, 

T a =l a ' , 2 a '=2 ao ' ; f*°' -f*° 


(13) 


From an examination of equations (ll) and. (3-2) , one can realize 
that the matrices that will he needed are (cp), (0), and (cp) ( 0) which are 
proposed "below in tabular form 



\ 

I a ’ 

—a 1 

2 

?■ 

\ 

- 


> 


I s 

1 

0 

0 


—a— a’ 
1 *1 

—a —a' 
1 -2 

—a —a’ 

1 -3 

(?) : 

2 s 

0 

ccp 

scp 

= 

2 a *T a ' 

—a —a* 

2 *2 

— a — a * 

2 *3 


— & 

3 

0 

-scp 

cep 


- 3 a .i a ‘ 

—a —a* 

3 *2 

—a — a ' 

3 *3 


V 


2*°’ 

3 a o’ 

\ 





I 8 -' 

C0 

0 

-S0 


I a ‘ -I 80 * 

l a ’ •2 a °' 

l a * .380* 

(0) : 

2 s ’ 

0 

1 

0 

= 

—a* — ao* 

z -1 

2 a ’ .a 80 ’ 

2 s ’ .3 s0 ’ 


t 

60 

0 

C0 


2® •q a ° 

3 a ’ . 2 ®°’ 

3 a ' -I 00 ' 


\ 

r ao T 

2 ao ' 

r* 

\ 




r a 

C0 

0 

— S0 


^a-ao’ 

l 8, • 2 80 ' 

l a -l 80 * 

(cp) ( 0) : 

—a 

2 

scps0 

ccp 

scpc0 

= 


a — ao 1 

2 • 2 

— a — ao* 

• 3 ° 


—a 

3 

ccps0 

-scp 

ccpc0 


- 3 a-a°’ 

3 & . 2®°* 

-a -ao' 
3*3 J 


( 1*0 


Now, returning to equation (ll) one can solve for q>, $, and by forming 
the appropriate dot products, then consulting the tables in (l4) : 


—a — a — a —a —a — a * — cl —a * ■ — ——a 

4 > = pi -1 + q 2 -1 + r 3 -1 -62 -1 -^3 -1 


( 15 a) 


= p+ijrs0 
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6 = pl a -2 a '+q2 a -2 aI +r3 a -2 a ' -<pl a -2 a ' -t3 & -2 a ' 
= qccp-rsq> 


( 15 b) 




t = pl a *3 ao ’+q.2 a >3 ao, +r3 a -3 ao, -cf)l a -3 ao -©2 a ' -S 80 ’ 

= -(p-9)s8+qscpc0+rccpce 
= fs 2 6+qscpc6+r«pc0 
f(l-s 2 0) = fc 2 0 = (qs<p+rcq))c0 
■fc0 = qscp+rccp 


>(15o) 


The term (p - qi) in equation (15c) was eliminated by using equation (15a). 
The angular velocity of the earth coordinates is 


fl 6 = I r +n| r 2 r +2s r 3 r = Q^'+yf 0 ' 


The primed coordinates here are the same sort of intermediate coordinates 
as those in equation (ll). The procedure just completed for &, i|r 
could he repeated here. But since^the only changes are the names of the 
components and the fact that <£=£>= 0, § and f c an he written down 

just hy comparison with equations ( 15 ): 


9 = p+fsS 
9 = qccp-rscp 
fed = qs<p+rccp 



$c8 = 


a 


ie 


n 


2e 


n e 

“se 


+ ise 



A final step remains in the geometry, the determination of fi 6 . This 
vector is obtained from the knowledge of the airplane' s velocity: 

v = R = — H-iifxR 

5t 
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The term R expresses the total rate of change of position of the air- 
plane, that is, its velocity. When this velocity is expressed in the 
rotating earth coordinates, however, it takes the explicit form shown in 
the expression on the right. The last term of this expression was derived 
in equation (7). It gives the full expression for the velocity in rotating 
coordinates if the vector R is constant in length. The other term in 

the expression above, — , accounts for changes in length of 5. Since 

St 


R = R1 


— — e — e — e — e Sr — e — e 

v = viel +v 2e 2 +v 3e 3 = 1 +Rftxl 

at 


or 


■= I s M +H ( %e 5 e -n 2 e3 e ) 

St 


v - SR 
Vle '^ 


v 2 e — R&3e = R^c0 
v 3e =“Rfiae = -R0 


These expressions complete the geometry. 

Although some of the transformations obtained along the way could 
have been obtained more easily without the use of the matrix method, this 
method is very valuable in providing a simple way of getting the complete 
expressions required in the problem. Furthermore, once a person gets the 
feel of the matrix method he appreciates its dual advantages of being a 
straightforward method of getting the correct expressions he wants and 
of providing a simple and clear view of the transformations required in 
the problem. The summary below, which actually was the analysis of the 
problem used before any computations were made, shows in simple fashion 
the various relationships involved in specifying the geometry: 
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[present airplane coordinates] 

= (cp)(e)('>|f)[ initial airplane coordinates] 

= ( 9 ) (sKtK -90°) [initial earth coordinates] 
= (a) (-{3} [present wind coordinates] 

[present earth coordinates] = (©)(¥) [initial earth coordinates] 


To express V = vl in earth coordinates: 

= (e)(T)[(9)(8Kt)(-<^)]- l (aH-p)^ 


The subscripts on the vectors indicate the coordinate system in which 
their components are to be expressed. To express g = (-g)! 6 in airplan e 
coordinates : 

la = |(®)(^)E(9)(0)(t)(-9O°)] -1 j- X g e 


To determine the rotation of the earth coordinates: 

v = I® p- + RO e xI e 

dt 


Position of a Bomber Relative to an Attacking Airplane 


The determination of the position of a bomber relative to an attacking 
airplane is a geometry problem that has to be set up in order to study 
weapon fire-Control systems by means of an analog computer. The final 
result sought in the geometry is the orientation of the line of sight and 
the range of the bomber from the attacker. As in the first example, with 
the proper precautions, a coordinate system can be assigned so as to make 
the line of sight correspond in direction to a unit vector of the coordi- 
nate system. With this assignment made, one seeks the rate of rotation 
of the line -of -sight coordinates required to make R = Rl^. Since the 
line-of -sight orientation relative to the airplane is desired, its rate 
of rotation will be expressed as the airplane rate plus a rotation rate 
relative to the airplane. Furthemore, since the study of automatic 
fire-control systems involves the use of a tracking radar, it will be con- 
venient to express the relative position of the line of sight by means of 
the same kind of angles by which the relative orientation of the tracking 
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line or heading of the radar antenna is expressed. The precaution 
mentioned above concerning the assignment of a coordinate system to the 
line of sight can be specified in a way to bring about the correspondence 
between the tracking line and line-of -sight systems. The condition will 
be: since the antenna does not bank relative to the airplane, neither 

will the line -of -sight coordinates. 

The velocity of the bomber, v* 3 = v' 3 !^ 3 , is assumed known. The air- 
plane velocity given in wind axes is V s = v a l vr £ angle of attack, a, 
sideslip, p, and the airplane angular velocity, 2 a = pl a +q2 a +r3 a , are 
assumed to be obtainable from the solution of the airplane equations of 
motion. One can obtain the angular velocity of the line of sight, f 2 2 , 
by finding its components from the equation 

v 2 = vb -va = I 2 X I 2 (17) 


As before, the subscripts indicate the coordinate system in which the 
vectors are to be expressed. Equation (17) implies that the bomber and 
airplane velocities are to be expressed in line-of -sight coordinates. 

The transformation from wind coordinates (in which v = v 1 ) to 
line-of -sight coordinates can be achieved by 


[line-of -sight coordinates] = (E) (A) [airplane coordinates] 

= (E) (A) (a) (-p) [wind coordinates] 



Sketch (h) 


( 18 ) 


where (A) performs a rotation about 
3 and (E) performs a rotation 
about 2 1 as shown in sketch (h). 

The nomenclature and order of 
rotation correspond to usage in 
some present fire-control systems. 

The magnitude of these angles E 
and A specifying the relative 
orientation of the line of sight 
is the final desideratum of the 
geometry. 

If the bomber velocity coordi- 
nates are first transformed into 
wind coordinates, the rotation matrix 
in equation (18) will complete the 
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t rans formation to line-of -sight coordinates. It is convenient to trans- 
form the ‘bomber velocity coordinates first into the initial attacker wind 
coordinates, then transform them into present wind coordinates. The 
transformations can follow the scheme: 


[bomber coordinates] = (©^(Y^ 3 ) [initial wind coordinates] 
[present wind coordinates] = (® v ) (&*) (Y* 0 [ initial wind coordinates] 


which gives the transformation 


[present wind coordinates] 

= ('S' w ')( 0 w )(l£ w )(-lt ) )(-Q^ l )[boiiiber coordinates] (19) 


Okie property of rotation matrices that (y ' 3 )"’ 1 = (-ijrk) has been used in 
equation (19). Now, since (Y w ) and (-Y* 3 ) are rotations about the same 
axis, (Y^X-Y^) can be written (Y) where the angle Y = Y”'- Y* 3 . 

The complete problem now has the form: 








— 


— | 



— - — 


■ 


v ll 


&t/dt 



vb 


v& 



= 

2l 

= (E)(A)(a)(-p)< 

(®)(e)(y)(-e'b) 

0 

- 

0 


Y 3l 


- Ra s\ 



0 


0 


— 


— — 



— 



y 


( 20 ) 


The superscripts in $ w and 0 W can be dropped without ambiguity. The 
quantities i, 0 , f™" are to be obtained from H a . The angular rates, 

§b, Y* 3 , specifying the bomber motions are assumed known; E and A are 
to be obtained from the relative angular velocity II = H -H 8 ". The 
matrix (a)(-p) is given in equation ( 10 ). The matrix (E)(A) is obtained 
from (a)(-p) by substituting E for a and A for -p. Then 
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~cEcA cEsA -sE~ 


eacp 

-casp 

-sa“ 



(E)(A)(a)(-p) = 

0 

3 

1 


sp 

cP 

0 




_sEcA sEsA cE 


_sacp 

-spsa 

ca 




cE cAcacP+cE sAs (3 - 

sEsac)3 

-cE oAcas p+cE s Acp+sEs p sa 

-cEcAsa-BEccT 

= 

-sAcacp+aAsp 



sAcasp+cAcp 

sAsa 


_sEcAcacp+sEsAsp+cEsacp 

-sEcAcasp+sEsAcP-cEspsa 

-sEcAsa+cEca 


The matrix ( 0 ) (f) can "be obtained from the matrix (E)(A) by the appropriate 
literal substitutions. Then 



“1 0 0 “ 


C0C? 008? -B8 - 


”ceb 0 sG b_ 

(®)(0)(?)(-0 b ) - 

0 e® s® 


-8? C? 0 


0 10 


_0 -B® c®_ 


800? 80S? C0_ 


0 

0 

■% 

ta 

LI 


rcec , ?ce b +sese b 
-c® s?c0b+s® ( B0c?c0 b -O0s0 b ) 
8®s'rc0 b 4-c® ( S0C?C0 b -c8s8 b ) 


008 ? 

O®cY+B®B0S? 

-s®c?+c®s0sY 


c0c?se b -sQcQ b “ 

-c®s¥s8 b +B® ( S0c?s0' b +cece" b ) 

B®B?B0bK3® ( S0C?S0b+C0C0b) 


The quantities <j>, §, fy* , and A remain to be determined; 0 13 

and Y * 3 are assumed given. The three angular rates <5, 0, »nrl t by 
means of which.the airplane wind coordinates are specified, are entirely 
analogous to <p, 8 and i{r in equations (ll) and ( 15 ). Instead of p, q, 
and r, however, the wind coordinate rates an< 2 . fi v 3W must be 

used. These components can be found as follows. Let be the angular 

velocity of the airplane with respect to the wind axes. Then 


!p- w - = S^ a 


- w i a +n s " v 2 a +n a " v 3 a 
2 a sa 
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This equation is again equivalent to equation (ll) for the condition 
q> = cp = 0 . The individual expressions can he obtained from those on the 
right in ( 16 ) after the appropriate literal substitutions 


Si 


,a-w _ 


1 a 


= (3sa 


St 


,a-w _ 


2 a 


= a 


S^~ v = pea, 
sa 


V S £t““W 

The wind axis angular velocity, ? i = 17 -U has, therefore, the 
components in airplane coordinates sirL a > ^ ~ q-o> j 

flga" r +3 cos a- . These components are transformed into a components in 
tne wind axis system through the transformation inverse to (a)(-p). With 
that transformation done, < 1 , 0 , and ¥ v are expressed in terms of 
p, q, r, d, and 

The quantities A and E are found in similar fashion, being analogous 
to and a,, respectively. The relative velocity of the line of sight, 
il, equals Ifi-I 7 a . The two components Sl^ and are found from equa- 

tion ( 20 ); 71 a is expressed in line-of-slght coordinates by means of the 
trans formation: 12^ = (E)(A)T7 a . Then A and E are given by 

• “1 

Slgi = E = + p sin A-q cos A 

9 

SI 3 I = A cos E = Sr 3 i' - sin E(p cos A+q sin A)-r cos E 


Rotation Matrix Method and Gimbaled Mechanisms 


The matrix method gives an extremely simple method of analyzing 
gimbaled mechanisms because the angles in the matrices can be chosen to 
correspond to the gimbal angles of the mechanism. 

Consider the antenna depicted in sketch (i). The 
3 axis is free to rotate in bearings mounted on a 
base and carries a frame which can be called an 
outer gimbal. Call -F 0 the angle made by the 
rotation of the axis relative to the base. The 
2 axis is free to rotate in bearings mounted on 
the outer gimbal, and carries the radar antenna 
with it. Call -r^ the angle made by the 2 axis 
in turning with respect to the outer gimbal. 

These physical angles, each measuring a rotation 

Sketch (i) 
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about an axis of the mechanism, are the same as the analytical angles, 
each of which describes a rotation about an analytic axis. The analysis’ 
angle A is measured by -r Q , and E by -ri . 

An orientation gyro is another gimbaled. mechanism. Sketch (i) repre- 
sents a gyro if the 2 axis, - instead of carrying an antenna, carries another 
ginibal that supports the gyro rotor, and if this rotor has its spin axis 
along axis 1. An orientation gyro has the property that its spin axis 
tends to maintain its orientation in space no matter what the motion of 
the base. Then the orientation of coordinates in the base is related to 
axes fixed in the inner gimbal by (ref. 5): 


[base coordinates] = (f 0 ) (fi)[ inner gimbal coordinates] 


Customarily in the analysis of airplane motion, the present airplane 
orientation is described relative to an initial orientation by the angles 
i)f(yaw), 0, and cp . Hence 


[present airplane coordinates] = (cp)(0)(ijr) [initial coordinates] 

= (cp}(0) [yawed coordinates] 


By fixing a gyro to an airplane appropriately. To and Ti can measure cp 
and 0 directly. The correct arrangement for the measurement of cp and 0 
requires the 3 axis aligned with the negative airplane 1 axis, 2 with the 
airplane 2 axis, and 1 with the airplane 3 axis under the null conditions 
cp = 0 = 0 and (sketch(i)) F 0 = = 0. Airplane yaw is not measured. 

This gyro arrangement corresponds to the standard vertical gyro arrangement. 


Transformation of Moment of Inertia 


Because they offer a good example of the use of rotation matrices, 
the equations of transformation of moment and products of inertia and of 
the force and moment derivatives of airplane dynamics will now be derived. 
The special cases involving a transformation from one body axis to another 
through a constant 0-like angle are well known (e.g., ref. 6). But since 
transformation equations, if their derivation is not understood, can be 
misused, it is desirable to derive the transformation in general, and 
illustrate the general method by a particular and common example. 

The fundamental dynamical law containing the moment of inertia is 
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where p is the angular momentum vector, oj is the angular velocity 
vector* and (j) is used here to represent the moment of inertia matrix 



JxX 

J xy 

Jxz 

(J) = 

Jxy 

J yy 

Jyz 


_Jxz 

Jyz 

Jzz_ 


Jxx 

J xy 



p dr 
etc 


where p is the mass density of the body at location x, y, z, and dr 
is an element of volume. 

Hie clue to determining the law of transformation of (j) is the 
recognition that equation (21) is a general dynamical expression having 
the same form no matter in which coordinates it may happen to be expressed 
in any particular example. Independently of the coordinates in which they 
are represented, angular momentum and angular velocity are connected by 
equation (2l): 


Po = (J 0 )“o 


Pn = (J)«n 


( 21 ) 


Now we have seen in transforming airplane and bomber velocity that the 
components of a vector transform in the same way as do coordinate systems. 
Thus, if r is a constant angle of rotation about some axis so that 


then 


[new coordinates] = (r)[old coordinates] 

Va. = (DPo > ^n = (r)«b or = (r) _1 u n 

in = (r)Po = (rKJoKr)- 1 ^ = (J n F n 


then 
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The transformation sought has the form 

(j n ) = (r)(J 0 )(r)- x 

If (r) is a transformation through a 6-like angle, then 



cr 0 -sf 


Jxx Jxy Jxz' 


""cr 0 sf 


(Jn) = 

010 


J xy J yy J yz 


010 



sT 0 cr 


Jxz Jyz Jzz 


-sr 0 cr 



j xx c2r+j zz s2 r -2j xz srcr 

Jxy°F -JyzSF 

( Jxx - Jz z ) cr+Jxz ( c£ r' --s 2 r ) 

= 

Jxycr -Jyzsr 



J yy 

Jxysr +Jyzcr 


i Jxx- Jzz ) sFcr+ J xz ( c 2 r -s 2 r ) 

JxySf+JyzCr 

Jxx s2 D +Jz z ° 2 r +2 Jxz cr _ 


( 22 ) 


Transformation of Force and Moment Derivatives of 
Airplane Equations of Motion 


To determine the law of transformation of the force and moment deriva- 
tives occurring in the equations of airplane motion, it is necessary, as 
with the moment of inertia, to find the derivatives in an equation 'whose 
form remains constant under the transformation. If the dependent variables 
are the six components of linear and angular velocity, the equations can 
be written in vector form as follows: 


A[E]+(n)[h] = (MiHvMMaHTn+CMsHs] 

ot 

~ [p]+(ft)[p] = (M 4 )[v]+(M 5 )[n]+(M 6 )[6] 


( 23 ) 


where h is the linear momentum, p is the angular momentum, v is the 
linear velocity, £1 the angular velocity of the airplane, and 



5G 


NACA TN 3968 


33 


(ft) - 

"0 -r q 
r 0 -p 

> (Mi) - 

'Hi 0 & 

o y v o 

S (Ma) - 

0 Xq 0 
Yp 0 Yr 

J (Ms) - 

0 Xse 0” 
0 0 Ygr 


21 P Q. 


_2u 0 


0 Zq_ 0 


0 0_ 


[»] - 

8 a 

8 e 

; [7] = 

V 

V 

S (M*) - 

"0 Lv d 
Hi 0 Mw 

i (Ms) - 

0 Lr 
0 Hi 0 

; (Mo) = 

L8a 0 Lsr 
0 tfee 0 




V 


_0 Nv 0 


Np 0 Nr 


jBa 0 N6r 


The transformation law of each of the (Mi) in equation (23) is identical 
to that of the moment of inertia, and is derived in the same way. 

Glauert’s formulas (ref. 6) are obtained using the 0-like rotation matrix 
(r) as in equations (22). For example, replacing (j) in equations (22) 
by (Mi) one obtains for the new X'u the expression: 


X’u - X u cos 2 r-(X w +Z w ) sin T cos r+Zw sin 2 r 


CONCLUDING REMARKS 


The present paper has attempted to show how convenient it is to use 
rotation matrices when one is setting up the geometric aspects of dynamical 
problems. A rudimentary knowledge of certain aspects of matrix algebra 
provides a tool entirely adequate for the solution of these geometric 
problems. The compact notation of matrix algebra permits the clear view 
required for their straightforward solution. The detailed computation 
of the expressions required in a particular problem becomes a matter of 
routine, and can easily be checked for errors . No algebraic dexterity 
is required to determine the parameters involved because they are obtained 
by the same direct methods used in the rest of the problem. Many of the 
parts used in one problem can be saved and used in other problems. 

Finally, because planning and computing become distinct tasks, it is a 
simple matter for one to devise and investigate various paths to the 
solution of the geometric problems without perfo rming any computations. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 4, 1957 
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APPENDIX 

A DESCRIPTION OF NOTATION 


Problems whose answers are currently being sought on simulators have 
caused a need for a uniform and consistent notation. The need flows from 
two sources: frequently more than one quantity of a sort is used in a 

single problem; and a single quantity is often required in more than one 
coordinate system. The person who is studying the problem is certainly 
familiar with the meanings of the symbols he employs. But if, as often 
happens, he employs a new letter for each new quantity, and still new 
letters for some components of these quantities, his audience must have 
constant recourse to a glossary of symbols. Some people, to avoid this 
pitfall, use a shorthand of subscripts by which they distinguish between 
quantities of the same sort. Unfortunately, this shorthand, which is not 
notation, quickly becomes too cumbersome. 

The notation proposed for use in the simulation of dynamical problems 
has three fundamental parts. The base letter of a quantity refers to 
quality; it indicates the sort of thing described. Superscripts are used 
to distinguish two quantities of the same quality. Subscripts are used 
to specify the coordinate system in which a quantity is represented, and 
to name the components of a quantity in a particular coordinate system. 

For example, the symbol v designates a linear velocity; v"b is a particu- 
lar linear velocity, say, bomber velocity; signifies the representation 
of bomber velocity in some particular reference system, for example, in 
attacker wind coordinates. 

When modification of a base letter by superscripts or_subscripts is 
not required for clarity, it should not be used. Use of D for the 
angular velocity of the line of sight relative to the attacker exemplifies 
the simplification possible. In this case, of all the angular velocities 
considered one need not have a superscript. It is convenient that the one 
angular velocity thus singled out be either the one most used, or the one 
that otherwise would have the most complicated superscript. There are 
many places where modification of a base letter is superfluous. 

Certain letters are favored over others for certain sorts of quanti- 
ties because of custom, or because of their mnemonic value. It seems to 
be desirable to reserve Greek letters for angles and angular velocities. 

The angles $>, 0, Y have been used to describe the orientation of a body 
with respect to a fixed coordinate system. Other letters describe ori- 
entation with respect to moving coordinates. Each of the letters $, 0 , 
if refers to a particular sort of angle, that is, one generated by a 
rotation about a 1, 2, or 3 direction. It seems preferable to capitalize 
these letters simply for typographical reasons If they are to bear super- 
scripts. The use of the small letters cp, Q, and i|r to indicate the 
orientation of an airplane body coordinate system, and p, q, and r to 
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represent airplane angular velocity components in body coordinates is 
honored, by custom. It seems undesirable to change customary symbols when 
their use is not confus ing and when they adequately fit the description 
needed. 

Another type of notation, a schematic representation of rotation 
operations, provides a convenient formulation of the geometry involved 
in a problem. It is derived from the fact that resolvers are the physical 
means of performing rotations in analog equipment. 

Induction resolvers (ref. 7 ) have the property of converting input 
quantities Ao and B 0 to output quantities An and Bn where _ 

An =_Aoeos r+B 0 sin T and B n = -Aosin r + B 0 cos I*. That is, if A = 1, 

B = 2, and r = Y, the matrix (Y) and the resolver below perform 


In' 


lo 

lo 

+ 

> In 

2n 

= w 

2o 

2o 


» 2n 

_3n. 


-3o_ 

3o 


3n 


equivalent operations. The presence of the + sign and the short circuit 
in the resolver symbol should be noted. The + sign indicates which out- 
put quantity edntains the +sin ¥ term. For an angle whose positive 
sense agrees with the positive sense of the coordinate system, the + sign 
follows the indicated short circuit in the cyclic order 1, 2, 3, 1, 2, 
.... Using this resolver symbol, one can write the matrix equation 


. (a)(-p)[(»)(e)(Y)(-*Tg-’?®] 


in the form 



V 


b 

lb 


0 

0 



(V b -V°) 

(V b -V°) 


(V b -V°) 


la 

2a 


3a 


The latter form is more explicit than the former. It conveys enough 
information that a person can use it rather than the various intermediate 
equations to aid him in wiring a problem on a s imul ator. 
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